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On the Casimir entropy between 'perfect crystals' 
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We give a re- interpretation of an 'entropy defect' in the electromagnetic Casimir 
effect. The electron gas in a perfect crystal is an electromagnetically disordered 
system whose entropy contains a finite Casimir-like contribution. The Nernst 
theorem (third law of thermodynamics) is not applicable. 
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1. Introduction 

It is well known that fluctuation interactions at nonzero temperature are 
entropic in character, a prominent example being the critical Casimir effect 
in liquid mixtures close to a continuous phase transition (see Ref.l for an 
overview). The electromagnetic Casimir interaction is also associated with 
an entropy that determines its limiting behaviour at high temperatures 
and/or large distances. 2 ' 3 The Casimir entropy for two material plates has 
recently attracted much interest also for low temperatures, as for certain 
situations a violation of the third law of thermodynamics (the Nernst heat 
theorem) has been claimed. 4,5 This has been used to argue in favor of a 
description where the DC conductivity of the metallic plates is ignored. 
Although the result of this theoretical prescription provides a better fit 
to recent experiments, 6 the situation is, however, not satisfactory from the 
physical point of view. In addition, a similar analysis for an experiment with 
laser-irradiated semiconductors 7 leaves open the meaning of the threshold 
value above which the DC conductivity should be included in the theory. 

Much has been said about spatially dispersive mirrors where the third 
law is verified, due to the anomalous skin effect, 8 and where a continu- 
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ous cross-over from a dielectric to a perfectly conducting response has been 
found. 9 We focus in these proceedings on a strictly local framework, mainly 
for simplicity, but also to show that this case is thermodynamically con- 
sistent as well. We shall see that, indeed, spatial dispersion plays only a 
small role in the range of wave vectors that are relevant for current Casimir 
experiments. 

We take up the interpretation of Ref.10 where the nonzero Casimir en- 
tropy found as T — > was associated to two oscillators coupled via a third 
one. Following this idea, we consider two half-spaces filled by an ideal elec- 
tron gas, separated by a distance L, and provide a direct calculation of the 
entropy per area S(L) = Hm S(L, T) in one of the two field polarizations. 
This calculation highlights the following point: 

The Casimir entropy S(L) < results from the coupling between two 
systems that are not in equilibrium as T — > 0. They are filled with a frozen 
magnetization and, in the local limit, have separately divergent (bulk and 
surface) entropies that characterize the disorder of the electromagnetic con- 
figuration. The Casimir entropy is the correction to additivity when the two 
bodies are close enough for their frozen currents to be mutually coupled 
by the quasi-static magnetic fields that 'leak' through their surfaces. The 
Nernst theorem is clearly not applicable for this disordered system. The 
situation is quite similar to the 'ideal conductor' (in distinction to a super- 
conductor, see Ref.ll) that does not reach thermodynamical equilibrium as 
it is cooled, because its random bulk currents freeze. (See Ref.12 where it 
is argued how special this ideal conductor case is.) 

2. Casimir entropy from frozen medium currents 
2.1. Motivation 

We have analyzed in a recent paper 13 the overdamped field modes to which 
the unusually large thermal corrections to the Casimir force between met- 
als can be attributed. Substantiating previous observations, ' 15 we have 
interpreted the characteristic frequency, £l = D/L 2 = {^aL 2 )^ 1 , in terms 
of a diffusion equation (diffusion coefficient D) satisfied by the magnetic 
field and electric currents in a medium with DC conductivity a. 

There is no contradiction between thermodynamics and fluctuation elec- 
trodynamics in this case. Fields and currents induced in the metal are 
clearly damped and lose energy into the phonon bath, say. In equilibrium, 
however, this is compensated by field and current fluctuations that are 
created by the bath. This concept can be traced back to the Einstein- 
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Langevin theory of Brownian motion 16 and is also the very essence of the 
fluctuation-dissipation theorem. 17,18 Moreover, for the quantum field the- 
ory, the quantum (or zero-point) fluctuations of the bath variables are an 
essential tool to establish at all times the commutation relations for the 
field operators. 19-23 In the field theory considered in Ref.13, we dealt with 
overdamped modes: if the wave equation were homogeneous, its eigenfre- 
quencies would be purely imaginary, similar to free Brownian particles. The 
quantum theory of Brownian motion 24 provides a consistent scheme for the 
quantum thermodynamics of this damped system. In this setting, nonzero 
entropies and even negative heat capacities find a quite natural explanation 
(see, e.g., Refs.25,26). 

In the particular case of an ideal electron gas (or 'perfect crystal'), 
the diffusion constant D = D{T) = 0{T 2 ). As the temperature drops to 
zero, the diffusion-dominated modes of the electromagnetic field do not 
reach a unique ground state, but remain in the classical regime <C T. 
This motivates the present calculation where the electromagnetic Casimir 
entropy between perfect crystals is re-derived within a classical model. 

The basic ingredient are transverse modes that extend throughout the 
bulk of the gas: static current waves interlocked with a magnetic field. 
The magnetic fields associated to the bulk currents leave one medium, by 
continuity, and cross the vacuum gap to the other medium in the form 
of (transverse) evanescent waves. This coupling between the two media 
changes slightly the wave vector of each current mode. Summing over all 
modes, we get a non-zero change in entropy that depends, quite naturally, 
on the separation L. It represents the distance-dependent change per unit 
area of the (much larger) entropy of the two frozen bulk systems. 

2.2. Lagrangian and conservation laws 

We start with the Lagrangian density 

£=^+en«.A-l(VxAf (1) 

where the field £ describes the displacement of a charged fluid element, A 
is the vector potential, n a constant background charge density and e a 
coupling constant with units of charge. The current density is j = en| [see 
Eq.(2) below], so that | represents a velocity field. The first term in Eq.(l) 
is thus the kinetic energy (density), the second one a bilinear coupling, and 
the third one the magnetic energy. Note that we neglect electric fields here. 
This is consistent if we make the assumption that V • £ = V • A = 0. The 
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first equality ensures that the medium displacement does not produce any 
charge density, the second one is the Coulomb gauge. The variation of the 
Lagrangian (1) with respect to A gives the Faraday equation 

eni - — V x (V x A) = (2) 
Mo 

In addition, the displacement field £ is a cyclic variable, hence we get a 
conserved momentum field 

-^j- = 0, tv = nm£ + enA (3) 

There are two ways to implement this conservation law physically. 27 

(i) In a London superconductor, the current density is tied at all times 
to the vector potential, with the momentum 7r being zero: 

Tie 2 

superconductor: j = en£ = A (4) 

m 

The Maxwell-Faraday equation (2) becomes 

(A" 2 - V 2 ) A = (5) 

where the Meifiner-London penetration depth A is given by the familiar 
expression A~ 2 = [i ne 2 /m — £l 2 /c 2 (Q is the plasma frequency). Eq.(5) 
does only allow for solutions that start at the surface and exponentially 
decay into the bulk on a length scale A (or shorter). Except for a surface 
layer of thickness ~ A, the interior of the medium remains free of magnetic 
field: the Meifiner-Ochscnfcld effect. 

From the London equation (4), we can also conclude that the Meifincr 
effect is maintained in time-dependent fields (at least with sufficiently slow 
variations; a detailed analysis clearly goes beyond the simple model con- 
sidered here). For a given frequency component u), the 'dielectric function' 
of the London superconductor can be read off from the polarization field 
associated to £: 

P = en£=^- = -^E (6) 

leading to the so-called plasma model e(u>) = 1 — SI 2 /io 2 . There is no vio- 
lation of causality here, if we read Eq.(4) as a retarded response function 
between the current density and the time integral of the electric field (i.e., 
the vector potential). 

The option (ii) that complies with the conservation law (3) corresponds 
to an ideal conductor: 

<9j dA 
ideal conductor: — = and — — = (7) 
at at 
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which means that currents, once created, are not damped and that the 
magnetic field is static. The value of the conserved momentum it is not 
restricted otherwise. The Faraday equation (2) then yields the field A in 
terms of its source j. Note that the magnetic field is in this case tied to the 
current density, similar to the scalar potential and the charge density in 
Coulomb-gauge electrodynamics. 28 Let us switch to reciprocal space with 
wavevector q: the vector potential A q created by the current is 

A q = ^ (8) 



so that the Lagrangian (1) becomes 



£(A 2 + q"W W 



q 

where V is the quantization volume. The conjugate momentum becomes 

w « = 7( 1 + w) jq (10) 

2.3. Normal modes and entropy 

The normal modes of the effective Lagrangian (9) can clearly be chosen 
as plane waves, labelled by wave vector q and polarization index /i. The 
associated Hamiltonian, expressed in terms of the canonical momentum 
field, is then 

The (classical) thermodynamics of this system is determined by summing 
the free energies of the normal modes over the quantum numbers q, fi. For 
one mode, we find by calculating the classical partition function {(i = 1/T) 

i^ = -Tlog /d^exp^e^l 2 ) =|log^ (12) 

where e q = (V/nm)(l + A~ 2 q~ 2 ) -1 determines the mode's energy. The 
entropy of this polarization mode is 



S ^ = -^T=—2 l0g ^ + 2 (13) 



where the equipartition term +1/2 comes from the (3 in the logarithm. 

When we sum this over all modes, the entropy becomes to leading order 
extensive in the volume V of the medium. The part that depends on the 
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surface is calculated in the usual way. Consider two media (total volume 
V) with parallel surfaces of area A facing each other at a distance L and 
write the entropy of the total system in the form 



where s is the (intensive) bulk entropy density, S sur { is the entropy per area 
of one (isolated) surface and S(L) the Casimir entropy per area. We can 
read the latter as the deviation from additivity in the system of two media: 
it thus describes how the disorder (or information content) of the two plates 
is changed by the coupling across the vacuum gap of thickness L. 

The physical mechanism for this coupling is the penetration of magnetic 
fields through the medium surface, as allowed for by the electromagnetic 
boundary conditions. In the vacuum between the media, the fields satisfy 
the Laplace equation V 2 A = 0: for a given wave vector k parallel to the 
surface, they 'propagate' perpendicular to the surface (along the z-axis, 
say) as evanescent waves ~ exp(±fcz) where k — |k|. For a single surface, 
only solutions that decay into the vacuum are permitted. In the gap < 
z < L between two surfaces, even and odd solutions coshfc(z — L/2) and 
sinh k(z — L/2) can be constructed. This is illustrated schematically in Fig.l. 



Fig. 1. Illustration of standing waves at the surface of an ideally conducting medium. 
We plot the component of the vector potential tangential to the surface. Thin line: 
isolated surface, thick line: mode between two surfaces with even parity. 

Both the surface entropy and the Casimir entropy can be calculated 
from the phase shifts of standing wave modes (see Refs.29,30 for details). 
For the surface entropy, 



S = Vs + 2AS suli + AS(L) 



(14) 




OO 




(15) 
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where q = (k, k z ), and the mode functions in the medium (z < 0) are 
proportional to e lk r " sin(fc z z + 9^) with = (a;, y) the coordinates parallel 
to the surface. From this form, we can also read off a 'reflection coefficient' 
for (time-independent) waves from within the medium, r p = — c~ 2l9f \ From 
a physical point of view, we can interpret the phase derivative in Eq.(15) as 
a density of modes in q-space, more precisely, its change due to the surface. 

For an isolated interface, the usual matching of the component of the 
vector potential tangential to the surface and its derivative at the interface 
yields in TE-polarization (current perpendicular to the plane of incidence 
spanned by k and the surface normal) 

rTE = fcrl' tan ^ E = -^ (i6) 

Here, k = |k| gives the decay constant of the evanescent wave on the vac- 
uum side. In the TM-polarization (current in the plane of incidence), the 
current has to satisfy the boundary condition lim j z {z) = to avoid the 

z^O 

build-up of a surface charge sheet. (That case would require electrical field 
energy in the Lagrangian (1) and is best described within a spatially disper- 
sive model. 29 ) This boundary condition immediately leads to r^M — 1, and 
there is no phase shift. The surface entropy in TM-polarization hence van- 
ishes, while it is logarithmically divergent at large q in the TE-polarization: 
•Ssurf ~ — (87tA 2 ) -1 log(g c A) with a short-range cutoff q c . One needs a non- 
local description of the material response (spatial dispersion) to get a finite 
result, see, e.g., Ref.31 for the surface self-energy. 

For the Casimir entropy, a local calculation is sufficient, as we shall see 
now: the reflection phases for even and odd modes in the vacuum gap are 
found as (we henceforward suppress the TE-polarization label) 

k kL k kL 

tan6U on (L) = — -^coth— , tan6» odd (L) = — ^ tanh — (17) 

The entropy per area for the two-surface system, 2S' sur f + S(L), is then 
given by Eq.(15) with 9 replaced by 6 even (L) + # dd(£)- Subtract twice the 
single-interface phase shift and calculate the quantity 

I _ r 2 Q -2kL 

cxp 2i[0 cvcn (L) + 9 odd (L) 29] = 1 _ {r , )2c _ 2kL (18) 



as can be checked with straightforward algebra. The Casimir entropy from 
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TE-polarized bulk currents becomes [combining Eqs.(13, 15, 18)] 

oo oo 

„ /r . fkdk fdk z / /3VA 2 q 2 \ d T , , 2 2feiN 

S ( L ) = / ~~z — / -r^ log —. , 9 „, - 1 Tr^-Im log (1 - r 2 e _2fei ) 

v ; J 2tt 7 2tt V 27rnm(l + A 2 q 2 ) / dk z 6 v ; 
o o 

(19) 

The terms independent of k z in the entropy per mode are irrelevant: after a 
partial integration, the integrated terms vanish because for k z — > 0, oo, the 
reflection coefficient r becomes real. Manifestly, short-wavelength modes 
with 2kL ^> 1 are suppressed, and a local theory is sufficient unless L 
becomes comparable to the length scales typical for spatial dispersion (mean 
free path, Debye-Hiickel screening length, Fermi wavelength). 

2.4. Calculation of the entropy 

Integrating Eq.(19) by parts, we have to evaluate the integral: 

j 7 dfc, fc,log(l-r 2 e- 2 ^) 

k J 2tiA 2 (fc 2 + fc 2 )(fc 2 + F + A- 2 ) 1 ; 

— oo 

where we recall that r is given by Eq.(16) above. We have extended the 
integration domain to — co < k z < +oo, using the property r(—k z ) = 
[r(k z )]*. Observe that r(k z ), as a function of complex k z , satisfies |r(fc z )| < 1 
in the upper half-plane, that the integrand vanishes at infinity, and evaluate 
the integral by closing the contour. There are simple poles at k z — ik and 
k z = \{k 2 + \^ 2 ) 1 / 2 . At the first pole, the reflection coefficient (16) vanishes, 
and we get from the second one: 

4 = ^log[l-r 2 l(M )e- 2 % r P'( fc .0)= (^A-^ + fc (21) 

As it happens, the reflection coefficient r p \(k,uj) for electromagnetic waves 
from a plasma half-space [dielectric function after Eq.(6)] appears here, 
evaluated at zero frequency and in the TE-polarization. If Eq.(21) is inte- 
grated over k, we get the 'entropy defect' calculated in the Lifshitz theory 
of the Casimir effect using the local dielectric function of a 'perfect crystal' 
[see, e.g., Eq.(20) of Rcf.5]: 

S(L) = J k -^h = j k -^ log [1 - r 2 l( fc, 0) c- 2 ^] (22) 
o o 

A switch to the integration variable to y = 2kL shows that 
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where the scaling function f(L/X), plotted in Fig. 2, is dimensionless and 
normalized to unity in the limit L ^> A. Indeed, in this regime, one may 
expand r p \(k, 0) in powers of k to get the asymptotic series [for higher terms, 
see Eq.(20) of Ref.5]: 

L > A : f(L/X) « 1 - + 12^ + 0{\/Lf (24) 



f(LU) 




5 10 15 20 25 30 



Fig. 2. Scaling function f(L/X) for the Casimir entropy of for the ideal electron gas, 
denned by Eqs.(22, 23). Dashed: Eq.(24). A = c/Q is the London-Mcifincr penetration 
depth (plasma wavelength). The same f(L/\) governs the 'entropy defect' for the elec- 
tromagnetic Casimir effect between perfect crystals (i.e., Lifshitz theory with the Drude 
dielectric function and scattering rate j(T) = 0(T 2 )), see Eq.(20) in Ref.5. 



3. Conclusions 

We have analyzed the Casimir entropy for the ideal electron gas, in particu- 
lar the contribution of electric currents frozen inside the bulk. This system 
shows (electromagnetic) disorder, and the third law of thermodynamic does 
not apply in its orthodox formulation. We have recovered the 'entropy de- 
fect' (negative Casimir entropy at zero temperature) reported in several 
places in the literature. Its thermodynamically consistent interpretation is 
the measure of the change in the disorder of the frozen currents due to their 
interaction through quasi-static magnetic fields. 
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